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1 13.  ABSTRACT  (Minimum  200  worch) 

Abstract— K  general  formalism  based  on  the  Green's  function 
method  is  given  for  multielectrode  semiconductor  lasers.  The 
effects  of  both  spatial  hole  burning  and  nonlinear  gain  are  in¬ 
cluded  in  this  formalism.  An  effective  nonlinear  gain  is  intro¬ 
duced  by  taMng  into  account  the  influence  of  the  laser  structure 
and  the  associated  distribution  of  the  mode  intensity  along  the 
cavity  length.  The  results  obtained  for  Fabry-Perot  and  dis¬ 
tributed  feedback  lasers  show  that  the  effective  nonlinear  gain 
could  be  considerably  enhanced.  Affected  by  the  laser  struc¬ 
ture,  the  nonlinear  gain  has  a  different  power  dependence  than 
expected  from  material  considerations  alone.  By  including  this 
effective  nonlinear  gain,  the  frequency  and  intensity  modula¬ 
tion  properties  of  multieiectrode  semiconductor  lasers  are  stud¬ 
ied.  A  general  linewidth  expression  is  given  which  includes  con¬ 


tributions  from  spontaneous  emission  and  carrier  shot  noise.  It 
is  found  that  the  effective  o-factor  affecting  the  linewidth  is  in 
general  different  from  its  counterpart  affecting  modulation  and 
injection  locking  properties  due  to  spatial  hole  burning  and 
nonlinear  gain.  For  lasers  with  uniform  intensity  distribution, 
the  effective  o-factor  affecting  the  linewidth  increases  or  re¬ 
mains  constant  with  increasing  output  power  dep—iii«g  on  the 
model  used  for  the  nonlinear  gain.  For  X/4  phase-shifted  dis¬ 
tributed  feedback  (DFB)  lasers,  the  effective  o-factor  affecting 
the  linewidth  is  slightly  larger  or  smaller  than  that  for  uniform 
lasers  depending  on  the  value  of  the  normalized  grating  cou¬ 
pling  coefficient.  The  linewidth  due  to  various  contributions  is 
calculated  for  both  uniform  intensity  distributed  lasers  and 
phase-shifted  DFB  lasers. 
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Abstract— A  general  formalism  based  on  the  Green’s  function 
method  is  given  for  multielectrode  semiconductor  lasers.  The 
effects  of  both  spatial  hole  burning  and  nonlinear  gain  are  in¬ 
cluded  in  this  formalism.  An  effective  nonlinear  gain  is  intro¬ 
duced  by  taking  into  account  the  influence  of  the  laser  structure 
and  the  associated  distribution  of  the  mode  intensity  along  the 
cavity  length.  The  results  obtained  for  Fabry-Perot  and  dis¬ 
tributed  feedback  lasers  show  that  the  effective  nonlinear  gain 
could  be  considerably  enhanced.  Affected  by  the  laser  struc¬ 
ture,  the  nonlinear  gain  has  a  different  power  dependence  than 
expected  from  material  considerations  alone.  By  including  this 
effective  nonlinear  gain,  the  frequency  and  intensity  modula¬ 
tion  properties  of  multielectrode  semiconductor  lasers  are  stud¬ 
ied.  A  general  linewidth  expression  is  given  which  includes  con¬ 
tributions  from  spontaneous  emission  and  carrier  shot  noise.  It 
is  found  that  the  effective  a-factor  affecting  the  linewidth  is  in 
general  different  from  its  counterpart  affecting  modulation  and 
injection  locking  properties  due  to  spatial  hole  burning  and 
nonlinear  gain.  For  lasers  with  uniform  intensity  distribution, 
the  effective  a-factor  affecting  the  linewidth  increases  nr  re¬ 
mains  constant  with  increasing  output  power  depending  on  the 
model  used  for  the  nonlinear  gain.  For  \/4  phase-shifted  dis¬ 
tributed  feedback  (DFB)  lasers,  the  effective  a-factor  affecting 
the  linewidth  is  slightly  larger  or  smaller  than  that  for  uniform 
lasers  depending  on  the  value  of  the  normalized  grating  cou¬ 
pling  coefficient.  The  linewidth  due  to  various  contributions  is 
calculated  for  both  uniform  intensity  distributed  lasers  and 
phase-shifted  DFB  lasers. 


I.  Introduction 

Multielectrode  semiconductor  lasers  have  attracted 
considerable  attention  in  recent  years  (1],  [2].  This 
is  largely  motivated  by  the  worldwide  development  of  ad¬ 
vanced  communication  systems,  in  which  optical  sources 
with  narrow  spectral  linewidth,  wavelength  tunability,  and 
flat  frequency  modulation  response  are  required.  Al- 
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though  the  theory  of  such  lasers  is  well  developed,  some 
phenomena  limiting  the  laser’s  ultimate  performance  have 
not  yet  been  completely  understood.  For  example,  line- 
width  rebroadening  occurring  at  high  output  powers  [1], 
[2]  is  contradictory  to  the  standard  laser  theory,  which 
predicts  a  monotone  linewidth  decrease  with  increasing 
output  power  [3].  Various  factors,  such  as  side  mode  in¬ 
fluence  [4],  I  //noise  at  high  output  powers  [5]— [7] ,  etc., 
have  been  involved  to  explain  this  phenomenon.  For 
highly  single  mode  semiconductor  lasers  driven  by  a 
noiseless  curr  it  source,  much  attention  has  been  paid  to 
longitudinal  inhomogeneities  and  the  nonlinear  gain. 

Longitudinal  inhomogeneities  include  two  principal  as¬ 
pects:  the  structural  inhomogeneity  and  the  functional  in¬ 
homogeneity.  The  structural  inhomogeneity  results  from 
the  multielectrode  nature  of  the  laser  structure,  in  which 
an  optical  grating  can  be  introduced  in  one  or  more  sec¬ 
tions  and  the  injection  current  in  each  section  is  separately 
controlled.  The  functional  inhomogeneity  originates  from 
the  phenomenon  of  spatial  hole  burning  [8],  [9],  occur¬ 
ring  even  in  a  homogeneous  cavity  structure.  Its  existence 
is  due  to  the  fact  that  the  intensity  distribution  is  not  uni¬ 
form  in  most  laser  structures,  especially  in  phase-shifted 
distributed  feedback  (DFB)  lasers.  This  nonuniformity  in¬ 
troduces  a  nonuniform  carrier  density  distribution,  which, 
in  turn,  affects  again  the  field  intensity  distribution,  and 
so  on. 

Longitudinal  inhomogeneities  can  be  easily  included  by 
using  the  Green’s  function  method  initially  proposed  by 
Henry  {10].  In  this  method,  the  solution  of  wave  equation 
corresponding  to  a  point  source  (spontaneous  emission 
event  in  the  laser  case)  is  called  the  Green's  function.  For 
distributed  excitation  sources,  the  general  solution  is  ob¬ 
tained  through  a  spatial  integration  of  the  Green’s  func¬ 
tion  weighted  by  the  correspondent  excitation  source.  Us¬ 
ing  this  method,  the  spontaneous  emission  rate  coupled  to 
the  lasing  mode  is  related  to  the  laser  structure  flO).  This 
method  has  been  generalized  to  include  complex  phase- 
amplitude  coupling  effects  in  DFB  and  composite-cavity 
laser  structures  [  1 1 J.  More  recently,  Tromborg  et  al.  have 
included  spatial  hole  burning  by  considering  rhA 
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skian,  appearing  in  the  Green’s  function  formulation,  as 
a  functional  of  carrier  density  distribution  in  multiclcc- 
trode  DFB  and  distributed  Bragg  reflector  (DBR)  lasers 
[  |2 J.  This  treatment  enables  them  to  take  into  account  any 
nonuniform  carrier  density  distribution  due  to  the  pres¬ 
ence  of  multiclcclrodcs  and/or  due  to  spatial  hole  burn¬ 
ing.  However,  the  results  show  that  the  spatial  hole  burn¬ 
ing  alone  can  not  explain  iinewidth  rebroadening. 

Another  phenomenon  that  should  be  considered  is  the 
nonlinear  gain,  which  should  include  the  material  nonlin¬ 
ear  gain  (spectral-hole  burning,  carrier  heating,  etc.)  [13]- 
[16]  and  the  nonlinear  gain  introduced  by  nonuniform  lat¬ 
eral  carrier  density  distribution  [17].  The  material  nonlin¬ 
ear  gain  due  to  the  spectral-hole  burning  has  been  studied 
by  using  a  density  matrix  formalism  [13],  [14],  By  study¬ 
ing  its  influence  on  Fabry-Perot  type  lasers,  it  was  con¬ 
cluded  that  the  material  nonlinear  gain  can  lead  to  line- 
width  rebroadening  [15].  Recent  results  obtained  by 
Olesen  et  al.  give  a  contrary  conclusion  for  the  same  type 
of  laser  structure  [18].  An  effort  was  made  by  Tromborg 
et  al.  to  include  the  material  nonlinear  gain  for  lasers  with 
a  nonuniform  intensity  distribution  using  the  Grem’s 
function  method  [12].  In  their  analysis,  the  Green’s  func¬ 
tion  method  is  directly  applied  even  for  a  nonlinear  di¬ 
electric  constant.  Conceptually,  this  is  not  permitted  as 
ihe  method  is  only  valid  for  a  linear  dielectric  constant. 
Consequently,  the  Wronskian,  resulting  from  solutions  of 
wave  equation  for  a  lincar  dielectric  constant,  can  never 
be  an  explicit  function  of  photon  density  distribution. 

The  purpose  of  this  paper  is  to  give  an  analysis  of  lasers 
bv  including  both  spatial  hole  burning  and  the  nonlinear 
gain.  Different  from  previous  approachs,  the  material 
nonlinear  gain  is  taken  into  account  through  a  perturba¬ 
tion  method.  This  enables  us  to  treat  the  nonlinear  gain 
within  the  validity  of  the  Green’s  function  method.  Using 
this  perturbation  method,  an  effective  nonlinear  gain  is 
introduced,  which  depends  also  on  the  mode  distribution 
along  the  cavity  length.  The  spatial  hole  burning  is  in¬ 
cluded  through  a  power  dependent  carrier  density  distri¬ 
bution. 

This  paper  is  organized  as  follows.  In  Section  11,  the 
Green's  function  method  is  generalized  to  multielectrode 
lasers.  In  Section  III,  the  effective  nonlinear  gain  is  dis¬ 
cussed.  In  Section  IV,  a  small-signal  analysis  is  per¬ 
formed,  leading  to  expressions  of  modulation  transfer 
functions  and  spectral  Iinewidth.  In  Section  V.  results  for 
lasers  with  uniform  intensity  distribution  arc  given.  Re¬ 
sults  fora  phase-shilled  DFB  laser  is  presented  in  Section 
VI.  Finally  a  conclusion  is  given  at  the  end  of  the  paper. 

II.  Genekalized  Rate  Equation 

The  starting  point  of  our  analysis  is  the  propagation 
equation  in  the  frequency  domain.  We  shall  concentrate 
our  attention  on  the  longitudinal  axis  alone,  although  the 
same  principle  of  analysis  could  be  used  for  transverse 
tuning  devices  such  as  tunable-twin-guide  DFB  lasers 
[19],  The  electrical  Held  in  the  laser  cavity  is  thus  gov¬ 


erned  by  the  inhomogeneous  scalar  Helmholtz  equation 

[10]: 

V?£„(z)  +  &JN,  S)Ejz)  =  Fjz)  (1) 

where  V?  is  the  Laplacian  operator  for  the  longitudinal 
coordinate  z.  =  u/c  is  the  wavenumber  in  vacuum. 
<M(V,  S)  is  the  dielectrical  constant  depending  on  the  fre¬ 
quency  u,  carrier  density  is  N,  and  photon  density  is  S. 
Finally,  Fu(z)  is  the  Langevin  force  describing  the  spon¬ 
taneous  emission.  Equation  (1)  is  valid  for  either  Fabry- 
Perot  or  DFB,  DBR,  lasers. 

The  complex  dielectrical  constant  is  written  as 

eu(N,  S)  =  [/i  +  j{g  -  ttt)/( 2*o)]2  (2) 

where  n  is  the  refractive  index,  g  is  the  gain,  and  aL  the 
internal  loss.  The  dielectrical  constant  can  be  split  into 
linear  part  and  nonlinear  part  eNL: 

eJN,  S)  =  et(u.  N)  +  eNL(«,  N,  S).  (3) 

The  nonlinear  part  is  intensity  dependent.  The  inclu¬ 
sion  of  the  nonlinear  dielectric  constant  makes  the  theo¬ 
retical  analysis  somewhat  complicated,  as  the  Green's 
function  method  is  only  valid  for  the  linear  case  [20].  To 
overcome  this  difficulty,  the  nonlinear  part  of  the  dielec¬ 
tric  constant  is  moved  to  the  right-hand  side  so  that  (1) 
takes  the  form: 

V;£w  +  kott£w  —  Fu  —  k^6NL  Eu.  (4) 

It  should  be  stressed  that  the  nonlinear  part  is  being 
treated  nonperturbatively  here.  The  general  solution  o'  the 
scalar  equation  is  obtained  by  using  the  Green’s  function 
formalism  and  is  given  by  [20]: 

EU(Z)  =  -  \  Gjz,  z')kfcNLEu  dz' 

J(L> 

+  \  Gjz.z')FJz')Jz'  (5) 

J(L) 


where  the  integration  is  performed  over  the  total  cavity 
length.  Cv(z,  z')  is  the  Green’s  function  given  by  [20]: 


GJz.  z') 


ZM>)Z.(z< ) 
RV,  A Kz)) 


(6) 


where  z>  =  max  (z,  z’)  and  z<  =  min  (z,  z’),  Z+(z). 
Z_(z)  are  two  independent  solutions  of  the  homogeneous 
equation,  satisfying  the  boundary  conditions  for  the  left 
or  right  facet  and  W(u,  N(z))  is  the  Wronskian  of  these 
solutions.  The  Wronskian  is  a  functional  of  the  frequency 
u  and  the  carrier  density  distribution  N(z).  It  is  not  ex¬ 
plicitly  dependent  on  the  coordinate  z.  Our  approach  is 
different  from  that  of  Tromborg  et  al.,  in  which  the  Wron¬ 
skian  is  considered  also  as  a  functional  of  photon  density 
distribution  [12].  This  is  not  permitted  within  the  validity 
of  the  Green’s  function  method.  In  our  approach,  the  pho¬ 
ton  density  distribution  does  not  appear  explicitly  in  the 
Wronskian,  although  it  affects  the  carrier  density  distri¬ 
bution  through  the  spatial  hole  burning  effect  [8],  [9]. 
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For  a  laser  system,  it  is  important  to  consider  the  laser 
oscillation  condition.  In  the  linear  case,  this  condition  is 
expressed  by  setting  the  Wmnskian  of  functions  Z,  (z)  and 
Z_(z)  equal  to  zero.  These  two  functions  Z+(z)  and  Z.(z) 
are  then  identical  to  one  solution  2^(z),  satisfying  the 
boundary  conditions  at  both  facets  and  representing  the 
longitudinal  distribution  of  the  electrical  field  in  the  laser 
cavity.  In  the  present  case  where  nonlinear  dielectric  con¬ 
stant  is  considered,  the  above  conclusions  are  only  a  first 
order  approximation.  Within  the  validity  of  this  approxi¬ 
mation,  we  obtain  from  (5)  and  (6): 

Zo(z')tSLEudz'+  f  Zo(z')Fu(z')dz'. 

Zq(z)  J(L)  J(L) 


This  newly  defined  nonlinear  gain  takes  into  account 
both  the  material  and  the  structural  dependences  and  will 
he  discussed  in  detail  in  Section  III. 

Some  assumptions  have  been  made  in  deriving  the 
above  rate  equation.  First,  the  nonlinear  gain  term  is. 
strictly  speaking,  a  convolution  between  the  effective 
nonlinear  gain  and  the  field  complex  amplitude.  This  is 
replaced  by  a  simple  multiplication.  This  assumption  is 
valid  as  long  as  the  field  complex  amplitude  variations  are 
much  slower  than  the  optical  frequency.  Second,  the  in¬ 
fluences  of  nonlinear  gain  on  field  distribution  have  been 
neglected.  Our  theory  can  be  considered  as  a  first-order 
perturbation  theory.  A  more  complete  theory  would  take 
into  account  these  higher-order  effects. 


(7) 

As  the  RHS  of  the  above  equation  is  independent  of  the 
cooruinate  z,  the  only  solution  of  the  electrical  field  in  the 
frequency  domain  Ea(z)  is  of  the  type: 

EJz)  =  0uZo(z).  (8) 

For  small  perturbations  due  to  the  nonlinear  gain  and 
the  driving  noise  sources  (spontaneous  emission,  carrier 
shot  noise,  etc.),  the  Wronskian  is  expanded  around  the 
linear  operating  point: 

bW  f  bW 

W(w,  N(z))  =  -(<*>-  u»o)  +  I  t —  A/V(z)  dz  (9) 
bo)  J(L)  bN 

where  «0  is  the  emission  frequency  at  this  point  and  AN 
is  the  deviation  of  the  carrier  density  distribution  from 
this  point.  By  using  the  properties  of  Fourier  transform 
the  following  rate  equation  for/30(r)  is  obtained  [10],  [1 1]: 

^  =  j j(S  -  «o)  -  j  [L)  WnAN  dz  +  ~  GNL( N,  S)j 

•  /So  W  +  F go  (0  (10) 

where  WN  =  (bW / bN) / (bfV / bu) ,  u  is  the  actual  lasing 
frequency,  different  from  u0  due  to  the  nonlinear  gain, 
0o (0  is  the  slowly  varying  complex  amplitude  given  by: 

M)  *  j  ^  expj(u  -  u)t  du.  (II) 

Fg o  is  the  Langevin  force  associated  with  the  complex  am¬ 
plitude  in  the  time  domain.  The  spontaneous  emission  rate 
R,p  is  given  by  [10]: 

\  Z*g0nQn,pZdz  \  dz 

4«2  J'L)  J(U 

^  =  IdW/bul1  '  °2) 

The  integration  in  (12)  is  restricted  to  laser  sections  with 
gain  only  since  passive  sections  with  different  bandgap  do 
not  contribute  directly  to  the  spontaneous  emission  rate. 
The  term  GNL  is  the  effective  nonlinear  gain  given  by: 

C«.W.  S)  -  -  [u  Zfcw*.  S)  dz.  (13) 


III.  Effective  Nonlinear  Gain 

For  Fabry-Pcrot.  DFB  and  DBR  lasers  with  a  refrac¬ 
tive  index  and  gain  varying  smoothly  in  the  cavity,  it  can 
be  shown  that  1 12]: 


bW 

Jn 

dW 

bu 


kotig,, ,(-<*„  +  j)Zl(z) 
2k{)  f  w 

-  \  nn  Z^z)  dz 


(14) 

(15) 


where  .j*,/  is  the  differential  gain  and  aH  is  the  linear  ma¬ 
terial  linewidth  enhancement  factor  [3].  Inserting  (15)  into 
(13)  together  with  eNL  =  jngNL/kn,  the  effective  nonlinear 
gain  can  be  finally  written  as  |2I|: 

ZUz)gNL(N,  S)  dz 

J(L) 

CN|.((V.  S)  =  vx - - -  (16) 

\  Z\(z)  dz 

J(L> 


where  ,i»NI  is  assumed  to  be  complex.  Its  real  part  repre¬ 
sents  the  nonlinear  gain  and  the  imaginary  pan  the  non¬ 
linear  refractive  index.  Thus  the  effective  nonlinear  gain 
is  a  spatial  average  of  the  material  local  nonlinear  gain 
weighted  by  the  squared  field  distribution  rather  than  by 
the  intensity.  This  newly  defined  nonlinear  gain  takes  into 
account  both  material  and  structural  dependences.  It  sim¬ 
plifies  to  the  material  gain  for  uniform  intensity  distribu¬ 
tions.  The  effective  nonlinear  gain  has  the  same  origin  as 
the  longitudinal  spontaneous  emission  enhancement  fac¬ 
tor  K.  |10|.  |ll]:  the  different  longitudinal  mode  distri¬ 
butions  { Zjz ).  m  -  0,  I,  •  •  •}.  forming  a  complete  set. 
are  not  orthogonal  in  the  Hermitian  product  sense,  due  to 
the  presence  of  optical  gain  in  the  cavity  [22]. 

In  fact,  when  the  mode  distribution  forming  a  complete 
orthogonal  set  is  used,  as  is  assumed  in  the  classical  laser 
theory,  the  effective  nonlinear  gain  is  given  by  [23): 

f  |Z..(<-)|WW-  S)dz 

Jll.l 


|Z,,(D|:  dz 


CN,  (/V.  S)  =  r. 


(17) 
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Several  consequences  of  the  effective  nonlinear  gain  can 
be  predicted  by  comparing  (16)  and  (17): 

i)  Different  laser  structures  can  give  rise  to  different 
values  and  forms  of  the  effective  nonlinear  gain  for  the 
same  material.  This  is  due  to  the  dependence  of  the  field 
distribution  Z^z)  on  the  cavity  structure. 

ii)  As  the  field  distribution  generally  includes  a  spa¬ 
tially  dependent  phase,  the  material  nonlinear  gain  can 
result  in  an  effective  nonlinear  index  and  vice  versa. 

iii)  As  the  intensity  distribution  is  not  uniform  and 
changes  with  the  output  power  due  to  spatial  hole  burn¬ 
ing,  the  effective  nonlinear  gain  will  have,  in  general,  a 
different  power  dependence  than  the  material  nonlinear 
gain. 

In  order  to  understand  the  implications  of  (16),  we  have 
to  consider  a  specific  functional  form  of  the  nonlinear  gain 
gNL.  However,  the  functional  form  depends  on  the  mech¬ 
anism  responsible  for  the  nonlinear  gain  (spectral-hole 
burning,  carrier  heating,  etc.).  In  many  cases  of  practical 
interest,  one  can  assume  that  gNL  decreases  linearly  with 
the  photon  density  (S(z))  as: 

Snl  =■  -8lS(z)/Ps  (18) 

where  P,  is  referred  to  as  the  saturation  photon  density 
and  gL  as  the  linear  gain.  This  expression  is  obviously 
valid  at  low  powers  such  that  S(z)  «  P,.  P,  is  a  material 
parameter  that  depends  on  details  of  carrier  relaxation 
within  the  conduction  band.  An  expression  of  P,  can  be 
found  in  [15]  for  the  case  in  which  spectral-hole  burning 
is  the  origin  of  nonlinear  gain.  Its  typical  value  is  in  the 
range  of  3-6  x  10IA  cm"'.  Equation  (18)  assumes  that 
gNL  is  real  such  that  the  refractive  index  is  power  inde¬ 
pendent.  In  the  case  of  lasers  operating  away  from  the 
gain  peak,  the  nonlinear  gain  is  accompanied  by  index 
changes  (15]  that  can  be  included  through  a  complex  gNL. 
In  this  paper,  such  changes  are  neglected  by  treating  gNL 
as  purely  real. 

By  relating  the  photon  density  to  the  mode  distribution 
through  S(z)  =  Po|2^(z)  |2,  with  P0  a  spatially  indepen¬ 
dent  parameter  proportional  to  the  output  power,  the  GNl 
can  be  written  as: 

[  Zl(z)\Z>(z)\- dz 

p  J(L> 

Gnl  -  -  ~  «V*l - - - •  (19) 

l  Z\{z)  dz 

JlL) 

It  is  assumed  for  the  moment  that  spatial  variation  of 
the  linear  gain  due  to  spatial  hole  burning  in  the  laser  cav¬ 
ity  can  be  neglected.  By  comparing  this  effective  nonlin¬ 
ear  gain  with  the  nonlinear  material  gain  G”lT  corre¬ 
sponding  to  the  average  photon  density  in  the  cavity: 

GT  =  “  vngL  (  I Uz)  |2  dz/ L  (20) 

r,  J(L) 

a  correction  factor  for  the  gain  can  be  introduced,  which 
relates  the  effective  nonlinear  gain  to  the  material  nonlin¬ 


ear  gain.  By  using  (19)  and  (20),  the  correction  factor  C 
is  given  by: 

f  Z20(z)\Z0(z)\idz 

C  _  Onl  _ _ ^<L>  _ 

(21) 

The  correction  factor  C  given  by  (21)  shows  how  the 
laser  structure  can  affect  the  value  of  the  material  param¬ 
eter  P,.  Note  that  C  is  generally  complex.  The  real  pan 
represents  the  change  in  the  material  nonlinear  gain,  the 
imaginary  part  represents  the  contribution  of  the  material 
nonlinear  gain  to  the  effective  refractive  index.  The  latter 
directly  affects  phase  variation  of  the  optical  field  in  the 
laser  cavity  and  thus  contributes  to  the  frequency  chirp 
and  the  spectral  linewidth. 

We  have  calculated  the  correction  factor  for  different 
types  of  laser  structures.  Fig.  1  shows  the  result  for  a  Fa- 
bry-Perot  laser  with  one  rixed  facet  reflectivity  of  30% 
(R2  -  0.3)  and  a  varying  reflectivity  /?,  of  the  other.  The 
real  part  of  C  is  slightly  larger  than  unity  for  small  values 
of  reflectivity  (R\  <  10~4).  Beyond  this  value,  the  cor¬ 
rection  factor  keeps  at  unity.  The  imaginary  part  of  C  de¬ 
creases  with  increasing  facet  reflectivity  J?,.  The  value  of 
the  imaginary  pan  is  quite  small  in  comparison  to  that  of 
the  real  part. 

The  correction  factor  C  is  plotted  as  a  function  of  the 
normalized  coupling  coefficient  in  Fig.  2  for  a  conven¬ 
tional  DFB  laser.  Both  facets  of  the  laser  are  assumed  to 
be  AR-coated  (R|  *  R2  =*  0).  The  real  part  of  the  correc¬ 
tion  factor  increases  from  0.87  for  kL  *  1 .0  to  1 .3  for  kL 
=  5.0.  The  imaginary  part  of  the  correction  factor  changes 
from  negative  values  to  positive  values  with  increasing 
normalized  coupling  coefficient.  It  is  not  surprising  that 
the  effective  nonlinear  gain  in  Fabry-Perot  and  conven¬ 
tional  DFB  lasers  is  not  very  different  from  the  material 
nonlinear  gain  (corresponding  to  the  average  photon  den¬ 
sity),  as  the  intensity  distribution  in  these  lasers  is  rather 
uniform. 

The  result  for  an  AR-coated  X/4  phase-shifted  DFB 
laser  is  shown  in  Fig.  3.  The  real  part  of  the  correction 
factor  becomes  larger  than  unity  for  kL  >  1.25  and 
reaches  2  for  kL  *  4.0.  The  imaginary  part  of  the  correc¬ 
tion  factor  changes  sign  at  kL  *  1 .25  and  becomes  neg¬ 
ligible  for  larger  values  of  kL.  This  is  due  to  the  fact  that 
for  kL  =  1 .25,  the  field  intensity  distribution  is  nearly 
uniform  inside  the  cavity  [8].  For  kL  larger  than  1.25,  the 
field  intensity  is  more  concentrated  at  the  center  of  the 
cavity.  Otherwise,  the  field  intensity  is  concentrated  near 
the  two  fuccts  (8). 

To  evaluate  its  output  power  dependence,  the  effective 
nonlinear  gain  is  calculated  by  using  (16)  for  an  AR- 
coated  X/4  phase-shifted  DFB  laser  with  kL  *  3.0.  In 
our  calculations,  the  field  distribution  and  the  linear  gain 
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Fig.  1.  The  real  and  imaginary  parts  of  the  correction  factor  C  as  a  func¬ 
tion  of  facet  reflectivity  R,  for  a  Fabry -Perot  laser.  The  reflectivity  of  the 
other  facet  ffj  is  assumed  to  be  30%. 


Fig.  2.  The  real  and  imaginary  parts  of  the  correction  factor  C  as  a  func¬ 
tion  of  the  normalized  coupling  coefficient  kL  for  a  conventional  DFB  laser 
with  AR-coated  facets. 


Fig.  4.  (a)  The  real  part  and  (b)  the  imaginary  part  of  the  effective  nonlin¬ 
ear  gain  G,g/vt  as  a  function  of  output  power  for  a  phase-shifted  DFB 
laser  with  nL  =  3.0. 

mation  result  as  (18).  The  real  and  imaginary  parts  of  the 
effective  nonlinear  gain  divided  by  vc(GN L/vK)  arc  plotted 
as  a  function  of  the  output  power  in  Fig.  4(a)  and  (b). 
respectively.  The  material  nonlinear  gain  corresponding 
to  the  average  photon  density  is  also  given  in  Fig.  4(a) 
for  comparison.  It  can  be  seen  that  at  low  output  powers 
(<  I  mW),  the  effective  nonlinear  gain  and  the  material 
nonlinear  gain  exhibit  nearly  the  same  power  dependence. 
The  ratio  of  their  values  is  close  to  1 .5.  the  value  of  the 
correction  factor  C  given  in  Fig.  3  for  kL  =  3.0.  When 
the  output  power  increases,  the  intensity  distribution  be¬ 
comes  more  uniform  due  to  spatial  hole  burning  [8],  [9]. 
Consequently,  the  real  part  of  the  effective  nonlinear  gain 
approaches  the  material  nonlinear  gain  at  an  output  power 
of  about  10  mW  and  continues  to  decrease  for  higher  out¬ 
put  powers.  At  the  same  time,  the  imaginary  part  of  the 
nonlinear  gain  tends  to  change  sign. 


Normltacd  coupling  coefficient 

Fig.  3.  The  real  and  imaginary  parts  of  the  correction  factor  C  as  a  func¬ 
tion  of  the  normalized  coupling  coefficient  kL  for  a  X /4  phase-shifted  DFB 
laser  with  AR-coated  facets. 


IV.  Small-Signal  Analysis 
To  study  the  dynamic  and  noise  properties,  one  has  to 
use  the  rate  equation  for  the  carrier  density  in  the  laser 
cavity: 


are  calculated  by  including  the  spatial  hole  burning.  The 
material  nonlinear  gain  is  assumed  to  have  the  form  [IS]: 

Snl  -  *l/>/1  +  2 S(z)/P,  -  gL.  (22) 

Compared  with  the  expression  given  in  [IS],  a  factor  2 
is  added  in  (22)  to  give  the  same  small  signal  approxi- 


dN(z ,  f) 
dt 


=  J(z,  t)  -  R(N(z.  t)) 


-  vKg(z,  r)S(z,  f)  +  Fv(z,  f)  (23) 

where  /(z,  /)  =  7(z,  t)/(ed ),  J{z.  t)  is  the  current  density. 
e  is  the  electron  charge,  d  is  the  thickness  of  the  active 
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layer.  S(z,  t)  is  the  photon  density.  R(N(z,  /))  is  the  non- 
stimulated  carrier  recombination  rate,  and  F„(z,  t )  is  the 
Langevin  force  representing  the  carrier  shot  noise.  The 
photon  density  is  related  to  the  held  intensity  distribution 
by: 

su,  t)  =  c0|zbwi2m  Co  -  ( j(u  \Uz)\z  dz/L} 

(24) 

where  C0  is  a  normalization  constant,  and  P(i)  is  the  av¬ 
erage  photon  density  in  the  cavity.  The  assumption  that 
the  spatial  and  temporal  dependences  of  the  photon  den¬ 
sity  can  be  separated  is  a  direct  consequence  of  the 
Green’s  function  method  [see  (8)].  This  is  a  good  approx¬ 
imation  for  modulation  frequencies  much  lower  than  the 
mode  separation  frequency  of  the  cavity. 

The  above  analysis  defines  the  basis  of  our  formalism. 
In  the  following,  it  will  be  used  to  discuss  the  dynamic 
and  noise  properties  in  multisection  lasers.  It  is  more 
practical  to  convert  the  held  complex  amplitude  rate  equa¬ 
tion  into  photon  density  P  and  phase  0  rate  equations  by 
using  0o(r)  =  -Jp(i)  exp  (y'0(f)): 

dP  f 

—  =  2  WmLMz)  dz  P(r)  +  Gnl,(N.  S)P(i)  +  />(/) 

(It  Jit.) 

(25a) 

d<b  _  f 

—  =  -  u  -  \  WNr  A  N(z)  dz 

dt  J(L) 

+  jCNU(\.S)  +  F*(l).  (25b) 

The  subscripts  r  and  /  represent  the  real  and  the  ima¬ 
ginary  parts,  respectively.  The  stationary  solution  is  ob¬ 
tained  by  setting  d/dt  =  0  and  eliminating  the  noise 
terms.  The  obtained  w,  /\  N(z),  Z^(z),  etc.  arc  power  de¬ 
pendent  stationary  values  and  distributions  because  of  the 
effects  of  spatial  hole  burning  and  nonlinear  gain.  Small 
deviations  from  the  stationary  point  are  obtained  by 
linearizing  the  rate  equations.  If  AP(Q),  AN(z.  SI),  Ao>(ft) 
represent  the  Fourier  transform  of  deviations,  the  small 
signal  solution  is  given  by: 


Au(Q) 


J(L) 


U\Mz)\2q(z)\1  dz  +  Gnla/2 


jQ  -  GNLW/>  +  2 P  f  U2A(z)\Uz)\z  dz 
Jtu 

•  Ff  +  F+  +  f  HAz,  Q)(AJ  +  F„)  dz 

J(L) 

(26c) 

where  i4(z),  r„,  Uu  U2  are  given  by: 

=  O,  (f^  +  I ) 


(27a) 


1  =  dR(N(z))\ 

t*  am 


+  V,  ^  S(z)  (27b) 


„  WNr  -  1/2  dG^/dN 
u'~  jQ  T  1/r. 

W*  +  1/2  dGHU/ZN 
Ul~  JQ  +  1/r, 

3Gnl  _  VgdgNL(N,  S)/dNZUz) 
9N 


\  Z](z)  dz 

Jtu 


(27c) 

(27d) 


The  local  modulation  transfer  functions  denoted  by  Hp(z, 
0)  and  Hf(z,  0)  are  given  by: 


Hp(z,  Q)  = 


2U-, 


jQ/P  +  2  f  U2A(z)  |  Zviz)  |2  dz  -  Ghu., 

"(L) 


(28a) 


Hf(z.  Q)  = 


~U\ (jQ/P  —  GNLpr)  +  GNL/»,  U2 


jQ/P  +  2  1  U2A(z) \Zq(z)  |2  dz  —  GNu»r 
Jtu 

(28b) 


A  A'(c.  fl) 


Aj(z,  Q)  +  Fn  -  /t(z)|Zo(z)l2  AP(Q) 
jQ  +  1/r* 


(26a) 


zlP(Q)  = 


_ Fn _ 

jQ  ~  Gfnj>,P  +  2 P  1  U2A(z)  |  Z>(z)  |  dz 
Jtu 


+ 


Hp(z ,  Q)(AJ  +  Fn)  dz 


(26b) 


3Gnl 

Gnlp  - 


~  c*v* 


Zo(z)|Z0(z)|23gNL(N,  S)/dSdz 


(28c) 


In  obtaining  the  above  solutions,  it  has  been  assumed  that 
the  held  distribution  Z>(z)  remains  unchanged  even  in  the 
presence  of  perturbation  due  to  the  noise  and/or  modula¬ 
tion. 

In  Section  II  (Fourier  transform  from  u  to  /),  it  has  been 
assumed  that  the  carrier  density  distribution  A.V(z)  is  a 


ISO 
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constant  in  the  time  domain.  This  is  a  reasonable  approx¬ 
imation,  as  the  temporal  variation  of  AN(z)  is  much  slower 
than  the  optical  frequency  u,,  (—  10"  s"1).  By  contrast, 
the  same  A N(z)  is  considered  time  dependent  in  this  Sec¬ 
tion  (Fourier  transform  from  /  to  0).  This  is  not  incon¬ 
sistent  with  the  first  pan  as  long  as  0  «  w0  (typically  8 
<  10"  s-'). 

A.  Dynamic  Response 

The  dynamic  responses  can  be  obtained  from  (26)-(28) 
by  eliminating  all  noise  terms.  For  the  case  of  uniformly 
injected  single  electrode  laser,  we  have: 

AP(8)  =  t  Hftz,  Q)  dz  (29a) 

J(U  eV 

Au(8)  =  l  Hf(z,  Q)  dz  (29b) 

J(U  ev 

In  obtaining  (29),  we  have  used  A J  —  A l/(eV),  with 
V  the  volume  of  active  layer.  A  convenient,  parasite-free 
measure  of  FM  is  the  chirp-to-modulated  power  ratio 
(C PR)  [17],  From  (28)  and  (29),  the  CPR  is  found  to  be: 

CPR  =  2 7KP  =  4xP  ijQ  ~  GfiU,r?  +  G^iP/^) 


The  transfer  functions  AP , ,  AP2,  •  •  •  usually  depend 
on  modulation  frequency  8.  Pure  frequency  modulation 
or  pure  inicnsity  modulation  have  been  obtained  using  this 
method  for  multielectrode  distributed- Bragg- reflector 
(DBR)  and  DFB  lasers  |24|,  125|. 

B.  Phase  Noise  and  Spectral  Linewidth 

The  Langevin  forces  representing  the  spontaneous 
emission  F;,( 8).  FJQ)  arc  delta-correlated  [10]: 

<F,<8)FJ.( 8’)>  =  2Dxx.6in  -  8'),  X,  X'  =  P.  <t>. 

(33a) 

The  nonzero  diffusion  coefficients  are  given  by  [10): 

2  =  RIP/(2D\  2D„,  =  2PRhJV  (33h) 

where  /  ( =  P  V)  is  the  total  photon  number  in  the  laser 
cavity.  The  carrier  shot  noise  is  assumed  to  be  delta-cor¬ 
related  in  time  and  in  space.  The  correlation  relations  arc 
found  to  be  [12]: 

<Fv(c.  Q)FUz.  8')>  =  2 Dm(z,  Q)6(z  -  ;')5(8  -  8') 

(34a ) 


(30a) 

where  an  effective  a-factor  affecting  the  modulation  prop¬ 
erties  is  defined: 

\  t/,(z,  8)  dz 

,  J(L) 

<**  =  ~  - - •  (30b) 

\  U2(z,  0)  dz 

J(L) 

The  CPR  expression  (30a)  has  the  same  form  as  the 
classical  one  [17],  but  with  the  effective  parameters.  It  is 
interesting  to  note  that  the  newly  defined  a-factor  (30b) 
is  modulation-frequency  dependent  for  lasers  with  non- 
uniform  intensity  distributions. 

In  the  case  of  a  two  section  laser,  the  excitation  source 
is  the  modulation  current  density  A J(z)  =  A/,  /s\ ,  for  val¬ 
ues  of  z  lying  within  Section  I  and  AJ(z)  =  AI2/  s2,  for  z 
within  Section  II,  where  r,  and  s2  are  the  cross  section  of 
areas.  The  carrier  density,  the  photon  density  and  the  fre¬ 
quency  deviation  can  be  obtained  from  (26): 

AP(0)  =  A,,  A/,  +  AP2A/2  (31a) 

Au(Q)  =  i4*|A/)  +  A+2AI2  (31b) 

where  the  transfer  functions  AP ,,  AP 2,  •  ■  •  can  be  ob¬ 
tained  directly  from  (26)-(28).  Using  (31a)  and  (31b). 
conditions  for  pure  frequency  and  pure  amplitude  modu¬ 
lation  are  given  by: 

~  =  -  for  pure  FM  (32a) 

A/i  Ap  | 

—■  =  —  for  pure  AM.  (32b) 

A/i  A$  | 


2£>,v,v(z.  n)  =  2[i/,g(z)nI;,S(z)  +  R(N(z))]/y 

(34b) 

<FP(Q)F„(Z,  a ’)>  =  2 Dpn(z.  8)  6  (II  -  8’)  (35a) 

2Dpn(z,  8)=  -2i >Kg(z)n,pS(z)/V  (35b) 

The  frequency  variation  driven  by  noise  at  zero  frequency 
(8  =  0)  is  obtained  from  (26c): 


Au  =  F0-  ae"TF,/(2P)  +  \  Hf(z ,  0 )FV(C.  0)  dz. 

J(L) 

(36) 

The  effective  phase-amplitude  coupling  factor  for  the 
linewidth  is  obtained  from  (26c)  and  written  as: 


_  noise 

<*err 


J 

J(Ll 

,y  '^N|J| 

Nr  2dN  j 

tpA(z)\Zo\2  dz 

(  1 

Jil  l  | 

+  13CH 

29N 

T#/f(z)|Zo|:  dz 

Zj  Cni./', 


^  GsU'r 


(37) 


Using  ( 14)— ( 1 6)  and  neglecting  the  nonlinear  refractive 
index,  the  effective  a-factor  for  the  linewidth  can  also  be 
written  as: 


«or'c  =  -Real(x)/Imag(x)  (38ai 
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(a) 


(b) 

Fig.  5  The  interaction  scheme  in  a  semiconductor  laser  due  to  la)  photon  density  tlucluation  caused  by  spontaneous  emission. 

(b)  current  modulation  or  light  injection. 


JtLi 


[(“<*«  +  jfl /til Jx, PkMz)  ~  jCi,dg/dS]\Zo\-Zfc)  dz 


X  = 


I  Zkz)  dz 

Jll.1 


(38b) 


By  comparing  (30b)  and  (37),  we  found  that  the  effec¬ 
tive  a-factor  for  the  iinewidth  is  different  from  its  coun¬ 
terpart  for  modulation  at  zero-modulation  frequency.  This 
is  explained  in  Fig.  5.  In  the  modulation  case,  current 
modulation  creates  carrier  density  variations,  which  give 
rise  to  a  simultaneous  change  of  photon  density  and  phase 
through  the  change  of  gain  and  the  index.  The  a™*  rep¬ 
resents  the  ratio  between  variations  of  photon  density  and 
phase.  In  the  Iinewidth  case,  the  Langevin  force  Fp  due 
to  spontaneous  emission  gives  rise  to  a  change  of  the  pho¬ 
ton  density,  which  leads  to  a  carrier  density  variation.  The 
latter  introduces  an  additional  phase  noise  through  the 
;hange  of  the  refractive  index.  The  a-factor  for  the  line- 
*'idth  relates  this  additional  phase  change  to  the  Langevin 
farce  Fv.  Thus  these  two  effective  u-faciors  have  different 
Jcfinitions.  More  importantly,  current  modulation  and 
.ponianeous  emission  affect  laser  dynamics  differently. 
:or  instance,  an  increase  in  electron  population  due  to 
urrent  modulation  leads  to  an  increase  in  the  output 
>ower.  In  contrast,  an  increase  in  the  output  power  due 
o  the  spontaneous  emission  results  in  a  decrease  of  elec¬ 
tion  population  to  conserve  the  laser-oscillation  condi- 
ion. 

The  Iinewidth  is  given  by  the  value  of  the  power  spec- 
rum  density  of  frequency  noise  at  zero  frequency  divided 
iy  2ir.  By  using  (33)-(36).  the  final  expression  is  given 

iv : 


t  «• 

•  1  |  H,{z.  0)  |:20nn(c,  0)  dz 

J<Lt 

n»M*c  i* 

-  \  2Rcal(W,<:..  0))2D,.N<:..  0)  dz.  <3‘J) 

2irP  Jn.i 


The  first  term  represents  the  contribution  of  phase  fluc¬ 
tuations  due  to  spontaneous  emission  and  that  of  photon 
density  fluctuations  via  the  well-known  phase-amplitude 
coupling.  The  second  term  is  the  contribution  due  to  the 
carrier  shot  noise.  The  last  term  represents  the  contribu¬ 
tion  of  cross  correlation  between  the  photon  density  and 
carrier  density.  The  carrier  shot  noise  induced  Iinewidth 
has  also  been  pointed  out  by  using  other  methods  [26], 
1271. 

The  above  discussion  completes  our  formalism  on  laser 
dynamics  and  noise.  To  get  some  insight  into  the  influ¬ 
ence  of  nonlinear  gain,  the  theory  is  applied  at  first  to 
lasers  with  uniform  intensity  distribution.  Results  for 
phusc-shittcd  DFB  lasers  arc  given  thereafter. 

V.  Application  to  Lasers  with  Uniform  Intensity 
Distribution 

For  Fabry-Perot  lasers  with  high  facet  reflectivities  (/?, 
=  R:  *  0.3)  and  specially  designed  DCB  lasers  [28],  the 
field  intensity  distribution  is  nearly  uniform  inside  the 
cavity.  In  this  case,  the  spatial  hole  burning  effect  can  be 
neglected.  The  effective  nonlinear  gain  becomes  simply 
the  material  nonlinear  gain.  By  using  Ws  =  cKgj(- aH  + 
j) / 2  [11],  the  rate  equations  in  this  case  arc  written  as: 

dP 

—  =  vKgdANP(t)  +  GNLr(N,  P)P{t)  +  Fp(t)  (40a) 

d<t>  _  I 

—  =  u)  -  wq  +  ~  aHvxgdAN(z) 

+  j  GmuIN.  P)  +  F.0).  (40b) 

The  modulation  and  noise  properties  have  been  widely 
studied  using  similar  or  the  same  equations  previously. 
We  are  essentially  interested  in  the  Iinewidth  and  the 
phase-amplitude  coupling  factor. 
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A.  Linewidth  Calculation 

Using  the  same  method  as  the  previous  section,  the 
linewidth  is  written  as: 

Af  =  jj^(l  +  (“*?“)')  +  ^  |W/(0)|22Dnn 

noutc 

-  ^  2Real(///(0))2DpN  (4 1 ) 

The  effective  phase-amplitude  coupling  factor  is  given 
by: 


none  _ 
“eff  “ 


1  dgNu  _  1  dgNu/df 

gd  dN  gjVtTKg  +  Pdg/dP 
1  dgNtr  ^  1  BgsLr/dP 

gddN  SdvtrK  8  +  Pdg/dP 


TABLE  I 

Lasi  *  Pa* ami  i  ms 


I’.IIUIIH  h  l> 

Symlmls 

Values 

(‘jvtly  lent 

/ 

400  |iin 

Thickness  »»l  active  layer 

- 

0  15 

Width  *d  active  lavci 

- 

2  0 

h  licet  ivc  retractive  imlc* 

n 

1  50 

Group  mtkx 

hk 

T  5<i 

Internal  loss 

a.mt 

40  cm  ' 

Linear  linewidth  enhancement  lac  lor 

4  0 

Bragg  wavelength 

l  -M)  ***** 

Carrier  lifetime 

r. 

2  ns 

DilTcrcnli.il  gain 

Km 

1.5  •  10  “  cm" 

Transparent  carrier  density 

/V„ 

1  0  10“  cm  ’ 

Derivative  ol  low 

do  /  dN 

-2.5  10  "  cm' 

Derivative  of  retractive  indes 

< In/dH 

-I  D  10  cm' 

Spontaneous  emission  factor 

"v 

2.0 

Salutation  photon  density 

P. 

5  0  10“  cm  ‘ 

where  gNL  *  Gnl/ vn  and: 


where  t,  is  the  carrier  lifetime  corresponding  to  the  non- 
stimulated  carrier  recombination.  It  is  assumed  that  the 
index  nonlinearities  are  neglected  in  what  follows.  Using 
the  nonlinear  gain  expression  in  the  low  output  regime 
(18)  and  (42),  (43),  the  effective  a-factor  for  the  line- 
width  for  P  -»  0  is: 


noise  aH 

<*'«"  ~  — ' 


'  -  P.  /P,. 


t\g,,Tr 


(22)  is  used.  The  ratio  between  at(r  and  aelf0  is  found  to 
be: 

(46) 

The  ratio  increases  with  increasing  photon  density,  as  is 
shown  in  Fig.  6  for  various  values  of  r.  The  material 
a-factor  |scc  (48b)]  normalized  by  ac(To  is  also  shown  in 
the  same  figure  for  comparison.  It  can  be  seen  that  the 
increase  of  the  material  a-factor  is  more  rapid  than  that 
of  the  effective  a-factor  for  the  linewidth. 

The  frequency  tuning  efficiency  is  obtained  by  inte¬ 
grating  the  local  tuning  efficiency  over  the  cavity  length: 


a 5T  _  (I  +  p)(  1  +  r) 
“ciro  vl  +  Ip  +  r 


where  Pe  denotes  the  interband  saturation  photon  density. 
Using  the  values  given  in  Table  I,  Pc  is  estimated  to  be  3 
x  1014  cm"3.  The  typical  intraband  saturation  photon 
density  is  Ps  »  3-6  x  1016  cm"3.  Thus  the  ratio  PJP , 
is  typically  of  the  order  0.01,  leading  to  a  correction  of 
the  a-fac*or  of  the  order  of  1  % . 

This  result  shows  that  the  effective  a-factor  for  line- 
width  differs  from  the  linear  material  parameter  aH  even 
for  zero  output  power.  Mathematically,  this  is  due  to  the 
fact  that  the  first  order  derivative  dgNL/dP  is  not  zero. 
This  result  seems  surprising  since  the  nonlinear  material 
gain  becomes  negligible  in  the  low  output  regime.  How¬ 
ever,  the  effective  a-factor  depends  not  only  on  the  mag¬ 
nitude  of  the  nonlinear  gain,  but  also  on  its  first  order 
derivative.  The  latter  does  not  vanish  even  in  the  low  out¬ 
put  regime. 

As  a  consequence,  the  effective  a-factor  at  high  output 
powers  should  be  compared  with  aeffn  rather  than  with  the 
linear  parameter  «#/.  First,  it  is  assumed  that  the  nonlinear 
gain  has  the  same  form  as  in  a  two-level  system: 


gAN  -  No) 

1  +  p 


p  =  P/Ps. 


(45) 


In  this  case,  the  ae(r  is  cxuctly  the  same  as  aeiro.  Sec¬ 
ond,  a  more  accurate  nonlinear  gain  expression  given  by 


By  using  (41).  (46),  and  (47),  the  three  contributions 
to  the  linewidth  are  plotted  in  Fig.  7  as  a  function  of  the 
inverse  of  the  normalized  photon  density.  It  can  be  seen 
that  the  contribution  from  spontaneous  emission  is  dom¬ 
inant.  This  contribution  decreases  with  increasing  photon 
density.  The  contributions  from  the  carrier  shot  noise  and 
the  cross  correlation  increase  with  output  power  because 
of  the  increase  of  the  noise  diffusion  coefficients  Dw  and 
Dkp-  However,  their  values  remain  many  orders  smaller 
than  that  due  to  spontaneous  emission. 

It  is  ini|Hiiiaiit  to  note  that  the  thermal  effects  are  not 
included  in  our  analysis.  In  practical  cases,  thermal  ef¬ 
fects  arc  dominant  mechanism  for  the  frequency  tuning  at 
low  frequency.  Another  important  point  is  that  the  drive 
current  noise  has  also  been  neglected  in  our  model.  How¬ 
ever.  the  drive  current  noise  is  converted  by  nonlinear 
gain,  especially  by  thermal  effects,  into  phase  noise  |29| 
The  recent  result  on  tunable  lasers  showed  a  very  impor- 
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Fig.  6.  The  normalized  effective  a-factor  for  the  linewidth  as  a  function 
of  the  normalized  photon  density  for  a  uniform  intensity  distribution  laser. 
The  parameter  is  the  ratio  of  the  mterband  saturation  power  to  the  inuaband 
saturation  power.  The  normalized  material  a-factor  is  also  shown  for  com¬ 
parison 


Fig.  7.  Calculated  linewidth  due  to  spontaneous  emission,  carrier  shot 
noise,  and  their  cross  correlation  as  a  function  of  the  inverse  of  the  nor¬ 
malized  photon  density  for  a  laser  with  uniform  intensity  distribution. 

tant  linewidth  increase  and  lineshape  change  due  to  the 
drive  current  noise  [30],  [31]. 


tsi 


Nee— Baad  phew  density  (PfFt) 

Fig.  t.  The  normalized  effective  a-factor  for  the  linewidth  as  a  function 
of  the  normalized  photon  density  for  a  phase-shifted  DFB  later.  The  pa¬ 
rameter  is  the  normalized  coupling  coefficient  nL.  The  normalized  a-factor 
for  the  linewidth  in  a  Fabry-Perot  laser  is  also  shown  for  comparison. 


This  is  exactly  the  material  linewidth  enhancement  fac¬ 
tor  in  nonlinear  gain  cases.  The  same  parameter  a£r*  ap¬ 
pears  in  the  technique  of  injection  locking.  For  every 
model  of  nonlinear  gain,  is  thus  always  larger  than 
the  linear  a-factor  aH.  The  measurement  results  recently 
made  by  Nakajima  et  al.  on  a  strained  quantum-well  DFB 
laser  confirmed  the  increase  of  this  parameter  with  in¬ 
creasing  output  power  [34]. 

It  is  thus  not  appropriate  to  measure  the  power  depen¬ 
dence  of  a for  the  linewidth  by  using  the  CPR  or  injec¬ 
tion  locking  method,  since  such  experiments  measure  the 
power  dependence  of  rather  than  that  of  a3r“.  Direct 

measurement  of  this  dependence  from  linewidth  is  also 
difficult,  as  many  other  factors,  in  addition  to  the  spon¬ 
taneous  emission,  contribute  at  the  same  time  to  the  li¬ 
newidth  [see  (39)]. 


B  Measurement  on  the  Power  Dependence  of  the 
Effective  a-Factor 

Various  methods  such  as  chirp-to-modulated-power-ra- 
tio  (CPR)  measurement  [32]  and  the  injection  locking 
technique  [33]  are  used  to  measure  the  effective  a-factor 
and  its  dependence  on  the  output  power.  It  was  usually 
recognized  that  the  a-factor  measured  from  those  methods 
is  the  same  as  that  for  the  linewidth.  Using  (30),  the  CPR 
for  lasers  with  uniform  intensity  distribution  is  given  by 
1171,  [32]: 


CPR 


Aw 
2  TAP 


_  iinhJ 

<*cil 

4t  P 


(48a) 


where  is  the  effective  a-factor  for  modulation  given 
by: 


a 


mod  _ 

ctf  ~ 


J_  <>Kni. 
dN 


(48b) 


VI.  Application  to  a  Phase- Shifted  DFB  Laser 

For  X/4  phase-shifted  DFB  lasers,  the  intensity  distri¬ 
bution  is  highly  nonuniform  due  to  the  presence  of  the 
phase  shift,  except  for  the  normalized  coupling  coefficient 
kL  in  the  neighborhood  of  1.25  [8],  [9].  The  field  distri¬ 
bution  is  obtained  by  solving  the  coupled  wave  equations 
and  the  carrier  conservation  equation  using  a  matrix 
method.  The  self-consistent  field  distribution  is  then  used 
to  calculate  the  effective  a-factor  and  the  linewidth. 

Fig.  8  gives  the  results  of  the  a-factor  for  the  linewidth 
normalized  by  aeffo  in  (44)  us  a  function  of  the  normalized 
photon  density  for  k L  —  1.0,  2.0,  and  4.0,  respectively. 
The  normalized  a-factor  for  the  linewidth  in  a  Fabry- 
Perot  laser  is  also  shown  for  comparison.  The  nonlinear 
gain  model  described  by  (22)  is  used.  It  can  be  seen  that 
for  all  values  of  the  normalized  coupling  coefficient  the 
effective  a-factor  increases  with  increasing  photon  den¬ 
sity.  This  factor  is  slightly. larger  than  the  a-factor  of  a 
uniform  laser  for  kL  >  1 .25  and  slightly  smaller  for  «/  < 
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Inverse  of  normaUad  photon  density  (Pi/?) 

Fig.  9.  Calculated  linewidth  due  to  spontaneous  emission,  earner  shot 
noiae.  and  their  cross  correlation  as  a  function  of  the  inverse  of  the  nor¬ 
malized  photon  density  for  a  phase-shifted  DFB  laser  with  *1  »  2.0. 

1.25.  If  nonlinear  gain  is  not  taken  into  account,  no  dif¬ 
ference  between  the  effective  and  the  linear  a-factor  is 
observed. 

Fig.  9  shows  the  three  contributions  to  the  linewidth  as 
a  function  of  the  inverse  of  the  normalized  photon  density 
for  a  phase-shifted  DFB  laser  with  kL  =  2.0.  The  contri¬ 
bution  due  to  spontaneous  emission  has  the  same  order  of 
magnitude  as  in  a  uniform  laser.  In  contrast,  contributions 
due  to  carrier  shot  noise  and  the  cross  correlation  are  much 
more  important  than  in  a  uniform  laser.  This  is  due  to  the 
fact  that  a  more  important  frequency  tuning  efficiency  is 
obtained  by  spatial  hole  burning.  However,  the  sponta¬ 
neous  emission  remains  the  dominant  contribution  to  the 
linewidth  in  the  usual  output  power  regime  (P/Ps  <  1). 

VII.  Conclusion 

In  this  paper  the  Green’s  function  method  has  been  ex¬ 
tended  to  include  spatial  hole  burning  and  the  nonlinear 
gain  in  multielectrode  semiconductor  lasers.  Different 
from  previous  approachs,  a  complex  effective  nonlinear 
gain  has  been  introduced  replacing  the  material  nonlinear 
gain.  Based  on  this  method,  a  generalized  rate  equation 
has  been  obtained.  Amplitude  and  frequency  modulation 
transfer  functions  have  been  given  for  single  and  multi- 
electrode  lasers.  A  linewidth  expression  has  been  ob¬ 
tained,  which  includes  contributions  from  spontaneous 
emission  and  carrier  shot  noise.  The  spontaneous  emis¬ 
sion  coupled  to  the  lasing  mode  and  the  effective  cr-factor 
for  the  linewidth  have  been  found  to  be  dependent  on  the 
stationary  carrier  density  distribution.  The  carrier  shot 
noise  leads  to  a  spectral  linewidth  proportional  to  the 
square  of  frequency  tuning  efficiency.  The  contribution  of 
the  correlation  between  photon  density  and  carrier  density 
depends  on  the  tuning  efficiency  and  the  effective  a-factor 
for  the  linewidth.  It  is  found  that  the  effective  a-factor 
affecting  the  linewidth  is  in  general  different  from  its 
counterpart  affecting  modulation  and  injection  locking 
properties  due  to  spatial  hole  burning  and  nonlinear  gain. 

For  lasers  with  uniform  intensity  distribution,  it  is  found 
that  the  effective  a-factor  for  the  linewidth  increases  or 
remains  constant  with  increasing  output  power  depending 


on  the  model  used  for  the  nonlinear  gain.  For  X/4  phase 
shifted  distributed  feedback  lasers,  the  effective  a-factor 
for  linewidth  is  slightly  larger  or  smaller  than  that  for  uni¬ 
form  lasers  depending  on  the  value  of  the  normalized  cou¬ 
pling  coefficient.  The  linewidth  due  to  various  contribu¬ 
tions  is  calculated  for  both  uniform  intensity  distributed 
lasers  and  phase-shifted  DFB  lasers 
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